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Abstract. Despite the growing evidence that long Gamma-Ray Bursts (GRBs) 
\ are associated with deaths of Wolf-Rayet stars, the evolutionary path of massive 

stars to GRBs and the exact nature of GRB progenitors remained poorly known. 
However, recent massive star evolutionary models indicate that — for sufficiently 
low metallicity — initially very rapidly rotating stars can satisfy the conditions 

■ for collapsar formation. Even though magnetic torques are included in these 
{SJ ' models, a strong core spin-down is avoided through quasi-chemically homoge- 
£S| , neous evolution induced by rotational mixing. Here, we explore for which initial 
CNl ■ mass and spin-range single stars of Z = Zq/20 are expected to produce GRBs. 
t-H ' We further find a dichotomy in the chemical structure of GRB progenitors, where 

lower initial masses end their lives with a massive helium envelope which still 

■ contains some amounts of hydrogen, while higher initial masses explode with 
C/O-dominated hydrogen-free atmospheres. 

^ ■ 

o : 

, 1. Introduction 

Long gamma-ray bursts (GRBs) are believed to o rigina te from rapidly rotating 
massive Wolf-Rayet stars (see Woosley &: Heger 2004 for a review). Interest- 
ingly, recent observations in dicate that GRBs occu r pref erentially in meta l poor 
environme nts (F ynbo et al. 120031: Co nselice et al. 2005; Gorosabel et al. 2005; 
Chen et al. l2005l : Starling et al. l2005ah . However, it has been questioned whether 
metal poor single stars are able to produce rapidly rotating WR stars as GRB 
progenitors, based on two reasons. First, stellar evolution models which in- 
clude magnetic torques indicate that the core loses too much angular momen- 
tum d uring the giant phase to pr oduce collapsar and GRBs (Heger, Woosley Sz 
Spruit 120051 : Petrovic et al. [2005) . Second, the lower the metallicity, the more 
difficult is the removal of the hydrogen envelope - without which jets from the 
central engine could not es cape f rom the star - even from very massive stars 
(see, however, Meynet et al. l2005h . 

Tw o rece nt independent studies by Woosley & Heger ((2005h and Yoon & 
Langer (2005) give a plausible solutio n to this problem. They considered so- 
called homogeneous evolution (Maeder ll987fl . where rotationally induced chem- 
ical mixing induces quasi-homogeneity of the chemical composition of the star 
throughout core hydrogen burning. In this case, single stars can become Wolf- 
Rayet stars without the need of stellar wind mass loss, and they avoid the giant 
phase that otherwise would cause a significant decrease of the core angular mo- 
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Fi gure 1. Evolution of 16 Mq stellar models at Z[ n it — 0.001, for differ- 
ent initial rotational velocities (Vmit/Vkepler = 0.3,0.4,0.5,0.8), in the HR 
diagram. 



mentum by magnetic torques. The above mentioned authors showed that such 
evolution can actually lead to the retention of enough angular momentum in the 
core to produce GRBs, if metallicity is sufficiently low (Z ^ Zq/W). Here we 
present stellar evolution models at Z = 0.001 which include rotation and mag- 
netic toques (Spruit f2002j). and systematically investigate in which conditions 
stellar evolution can lead to GRBs, via such an evolutionary path. 



2. Methods 



The stellar models are calculated with a hydrodynamic stellar evolution code, 
which includes the effect of rotation on the stellar structure, rotationally induced 
chemical mixing, an d the transport of angular momentum by magnetic torques 
(see P etrovic et al. 120051 and references therein). We follow Kudritzki et al. 
( 1989) for stellar wind mass loss of hydrogen ric h star s. Wolf-Rayet wind mass 
loss rates are adopted following Hamann et al. (.1995), but reduced by a facto r 
of 10, which corresponds to the recent estimates by Vink & de Koter (J2005). 
The effect of the enrichment of CNO elements at the stellar surface on the WR 
wind mass loss rate is also considered such that, with a given surface condition, 
WC stars with Xcno = 0.5 and Xh c = 0.5 have about 10 times higher mass loss 
rates than WN stars. A metallicity dependence of M oc Z - 69 and M oc Z 0,S6 
is a dopte d for hydrogen rich stars and WR stars, respectively, following Vink et 
al. (I2001T 1 and Vink & de Koter (l2005h . We do not consider overshooting in the 
convective region, but employ r ather fast semi-convection with an efficiency pa- 
rameter asEM = 1-0 (see Langer Il99lj ). Uncertainties due to the se ass umptions, 
and their effects on the results are discussed in Yoon &; Langer (200l). 
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Figure 2. Evolution of the internal structure of the model sequences of 
Mi„it = 16M Q and Z = 0.001, with Vinit/VK = 0.3 (left panel) and V in it/V K = 
0.5 (right panel). Convective regions are hatched. Gray shading denotes 
nuclear energy generation. 



Evolution of massive stars at Z = 0.001 



As discussed by Maeder (1983), the evolution of rotating stars can bifurcate 
according to the initial spin rate. While non-rotating, or slowly rotating stars 
evolve red wards, rapidly rotating stars evolve bluewards due to very efficient 
rotationally induced chemical mixing. Figure ^ illustrates this bifurcation in the 
H-R diagram, with 16 Mq models. Initially slowly rotating stars (Vinit/Vx = 
0.3 &: 0.4) transform into red giants with a massive extended hydrogen envelope 
(Fig. I2J). The resulting CO core mass is about 2.8 Mq, and the star is expected 
to explode as a Type II supernova leaving a neutron star as a remnant. On 
the other hand, the initially rapidly rotating stars (Mnit/Vk = 0.5 & 0.8), which 
become WR stars on the main sequence due to rotationally induced mixing, 
have only a tiny compact hydrogen envelope during the core He burning phase 
(Fig. 12). Importantly, such a tiny envelope cannot spin down the rapidly ro tating 
core, contrary to the case where stars evolve into red giants (Heger et al. 120051 : 
Petrovic et al. I2005D . Rather efficient angular momentum transport from the 
CO core occurs after the core He exhaustion, but the central cor e can retain 
enough a ngular momentum to produce collapsar (Woosley & Heger 120051 : Yoon 
&; Langer ,2005:). The CO core mass at the final stage is about 10.4 Mq, which 
is large enough to form a black hole. 

Not all stars which undergo quasi-homogeneous evolution end their life with 
a GRB. If the metallicity is too high (Z ;> Z q/10), mass loss result s in a too 
strong spin-down of the star (Woosley & Heger l2005l : Yoon & Langer |2005T > . At 
a given metallicity, homogeneously evolving stars of relatively low initial mass 
(~ 10 Mq) have a rather long CO core contraction time, and a rather massive He 
envelope after the He core exhaustion. In these stars, the CO core loses too much 
angular momentum by magnetic torques after core helium burning to produce a 
GRB. This imposes a lower initial mass limit (M m i n ) for GRB formation. We find 
M m j n ~ 20 Mm with slow semi-convection («sem = 0.01) as discussed in Yoon 
&: Langer ( 20051 ). and M m - m ~ 12 M q with fast semi-convection (cjsem = 1-0), 
in agreement with Woosley & Heger ( 20051 ) . 
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Figure 3. Final fate of rotating massive stars at Z = 0.001, in the plane 
of initial mass and initial fraction of the Keplerian value of the equatorial 
rotational velocity. The solid line divides the plane into two parts, where 
stars evolve quasi-chemically homogeneous above the line, while they evolve 
into the classical core-envelope structure below the line. Between the dashed 
lines is the region of quasi-homogeneous evolution where the core mass, core 
spin and stellar radius are compatible with the collapsar model for GRB 
production, while to both sides of it black holes are formed but the core 
spin is insufficient. This GRB production region is divided into two parts, 
where GRB progenitors do or do not possess a thick helium envelope (i.e. 
AMhc > 2.0M Q ). The dashed-dotted line in the region of non-homogeneous 
evolution separates Type II supernovae (SN II; left) and black hole (BH; right) 
formation, where the minimum mass for BH formation is simply assumed to 
be 30 M0 (see, however, Heger et al. 120031 for a comprehensive discussion on 
the issue). From Yoon & Langer l|200fift . 



The upper initial mass limit for GRB formation is imposed by different 
factors for different metallicities. At Z = 0.001, stars with Mj n it ^ 40 Mq expe- 
rience significant braking by rather strong mass loss, and cannot retain enough 
angular momentum in the core. Therefore, only those stars with 12 M Q <J 
-^init ^ 40 M are likely to produce GRBs, at Z = 0.001. This upper mass 
limit will decrease with increasing metallicity, due to stronger stellar wind mass 
loss. On the other hand, if the metallicity is very low (Z « 10~ 5 ), the braking 
induced by mass loss is less important even for very massive stars. However, 
homogeneously evolving stars with Mj n i t ^ 60 M yields CO cores more massive 
than 40 Mq and likely under go the pair-instability, preventing the formation of 
GRBs (Yoon & Lanarer 1200.^ . 
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Figure 4. Chemical structure of two GRB progenitor models, as examples 
for the two different expected progenitor types. The left plot shows a "WN- 
type progenitor" (M init = 16 M Q , Mnit/Vk = 0.5), with a final total mass of 
13.07 M , a CO-core mass of 10.38 M s , and hydrogen in its outer 0.2 M . The 
right plot shows a "WC-type progenitor" (M init = 30 M Q , V init /V K = 0.4), 
which ends up with 20.08 M and a CO-core mass of 16.88 M Q . 

In Fig. |3J we summarize the above discussion on the final fate of massive 
stars at Z = 0.001, in the mass-rotation plane, based on a grid of stellar models 
with di fferen t initial masses and spin rates (about 40 model sequences; Yoon & 
Langer l2006f ) . Here, fast semi-convection (osem = 1.0) is adopted. 

4. Properties of GRB progenitors 

Our models at Z = 0.001 predict two different types of GRB progenitors. Ini- 
tially less massive stars give less massive CO cores and more massive He en- 
velopes. As a result, stars with 12 M Q ^ M^t 25 M Q end their life as WN 
stars with rather massive helium-rich envelopes (AMh c <; 2.0 Mq). Those might 
be characterized as Type lb supernovae (Mazzali 2005, private communication; 
See Fig. 0). Remarkably, these stars also have hydrogen in their envelope (left 
panel in Fig.0J, which might be relevant to the high velo city HI absorption line 
observed in the afterglow of GRB 021004 (Starling et al. l2005bT ). On the other 
hand, more massive stars become WC stars in the end, which will explode as 
Type Ic supernovae (right panel of Fig. |3J). Our models also predict that GRB 
progenitors of "WC-type" are more compact, more massive, and have envelopes 
which are more enriched with a-elements, than those of "WN-type" . 

5. Discussion 

The initial spin rate distribution of massive low-metallicity star is unknown, 
and thus our models can not readily predict a GRB formation rate at Z = 
0.001. However, Langer & Norman ((20061 ; see also the General Discussion after 
Session G, in this volume) argue that if the majority of GRBs were restricted 
to metallicities below Z = Zq/W, about 5 percent of all massive stars with 
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such low metallicity should produce a GRB in order to obtain GRBs at a rate 
observed by BATSE. This requires that, at this and lower metallicity, more 
stars produce GRBs than stars are pre dicted to die as WR star due to stellar 
wind mass loss (cf. Meynet & Maeder [2005) . Consequently, a significant low- 
metallicity bias in GRBs would — if confirmed — not only be consistent with 
the quasi-chemically homogeneous evolution scenario of GRB progenitors. It 
would require that indeed the evolution of low metallicity massive star differs 
significantly from that of massive stars in our Galaxy, in support of the scenario 
outlined above. 
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